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Abstract 

Electron plasmas confined by an external magnetic field exhibit variations in a two-dimensional 
plane orthogonal to the confining magnetic field. A nonlinear fluid simulation code to investigate 
the properties of 2-D electron plasma wave turbulence in a nonuniform magnetoplasma has been 
developed. It is found that the presence of the density gradient convection by mean electric 
fields considerably influence the characteristic nonlinear interaction processes, such as the energy 
cascades and the cross-field electron transport. The initial random turbulent state evolves towards 
an intermittent state where forward cascade of vorticity coexists with an inverse cascade of electric 
potential fluctuations. The latter lead to the formation of large scale entities in 2D electron 
plasmas and can be alternatively understood by seeking exact nonlinear coherent vortex solutions 
in the form of a dipolar-like configuration. The energy cascades are governed typically by the 
Kolmogorov-like A; -5 / 3 spectrum. In agreement with the experimental observations, we find that 
the electron transport is improved significantly by the application of an externally imposed electric 
field. 



*Electronic address: dastgeer@u cr.edu 
^Electronic address: ps@tp4.rub.de 



1 



I. INTRODUCTION 



Electron plasma columns are routinely confined in experimental devices by means of 
an external magnetic field for hours. While such non-neutral plasmas [l| exhibit excellent 
stability property under ideal conditions, there exist critical issues that can potentially 
destabilize the electron plasma columns. For instance, the asymmetry of confining field and 
the presence of neutral gas particles are reported to often cause damages to the stability 
of the electron plasma columns Q]. During the process of rotation around the ambient 
field, the electron plasma particles collide with each other resulting thus in a cross-field 
electron transport. The cross field transport is therefore dominated by the electron-electron 
collisions in an electron plasma column, whereas collisions with the background neutrals do 
not contribute to the transport. However the advent of a rotating wall electric field imposed 
externally on the electron plasma column leads to an improvement in the stability |2j. An 
application of the externally applied radial electric field E in the presence of the confining 
orthogonal magnetic field B imparts a net poloidal drift vi « (c/i? 2 )E x B on the electron 
plasma column. Consequently, the plasma column rotates around the external magnetic 
field which supposedly enhances its stability. 

There exist morphological correspondence between the guiding center flow of electrons in 
a strongly magnetized electron plasma and the two-dimensional (2D) incompressible Euler 
fluid In a uniform plasma, the electron fluid dynamics is governed by the Navier-Stokes 
equation (NSE), which basically describe the evolution of the electron vorticity. The NSE 
equation admits interesting electron fluid turbulence behavior and localized vortex motions 

n 

in the form of clumps (positive vorticity) and holes (negative vorticity) [4J. The dynamics 
of electron plasma vortices in background vorticity distribution has been investigated both 

n nn 

analytically and numerically [4] as well as experimentally [5, 6]. 

In a nonuniform electron magnetoplasma containing equilibrium electron density gradient 
and fixed ion background, the dynamics of 2D electron drift modes (EDMs) is governed 
by the modified NSE [7]. Possible stationary solutions of the latter can be cast in the 
form of a global vortex pattern [7j, [8| and a double vortex Vortex structures 
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clouds are confined by dc electric and magnetic 
can parametrically excite electron zonal flows 
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ds. Furthermore, large amplitude EDMs 
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confinement in nonuniform magnetized plasmas. 

In this paper, we present the turbulent properties of nonlinearly interacting low-frequency 
(in comparison with the electron gyrofrequency) electron plasma waves (or EDMs) in a 
nonuniform magnetoplasma containing the electron density inhomogeneity. The manuscript 
is organized in the following fashion. In Sec. II, we discuss our 2D model equation for 
nonlinear EDMs and the possibility of frequency condensation due to the mode coupling 
process. Simulation results for a random initial turbulence state are described in Sec. III. In 
Sec. IV, we seek to understand the emergence of large-scale coherent structures by virtue of 
exact analytic solutions. The evolution of turbulent energy and the corresponding transport 
are discussed in Sees. VI and VII, respectively. The conclusions are contained in Sec. VIII. 

II. BASIC EQUATIONS 

Let us consider a strongly magnetized electron plasma in the presence of the equilibrium 
density gradient (drio/dx) across zB , where z is the unit vector along the z direction in 
a Cartesian coordinate system and B is the strength of the external magnetic field. The 
ions are supposed to form an immobile uniform neutralizing background. The perpendicular 
component of the electron fluid velocity in the presence of nonlinearly coupled low-frequency 
(in comparison with the electron gyrofrequency u ce = eB /m e c, where e is the magnitude of 
the electron charge, m e is the electron mass and c is the speed of light in vacuum) electric 
field E = — V0, where is the electrostatic potential, is 

^«^xV0+^-(| + , e -^Vi)v^, (1) 

where d/dt = (d/dt) + (c/B )z x V0, fi e is the electron collision frequency, and fi e is the 
electron kinematic gyro-thermal- viscosity. The nonlinearity in (1) comes from the nonlinear 
electron polarization drift. The electron advection term v ez d/dz does not appear in d/dt 
because of our 2D approximation. 

Inserting (1) into the electron continuity equation and using the Poisson equation, we 
obtain the modified NSE 

s + $o«-*vi) 

where uj pe = (4nn e 2 /m e ) l l 2 is the electron plasma frequency, ujh = (u 2 pe + ^ e ) 1//2 is the 
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upper-hybrid resonance frequency, and K n = —dno/dx. In the absence of nonlinear inter- 
actions, Eq. (2) yields 



u pe uj ce k y K n .a; 
u = k v u o + — -2T2 z "i- r ' ( 3 ) 



which is the frequency of the damped Doppler-shifted EDMs in the local approximation 
(viz. the wavelength being much smaller than the density gradient scalelength K~ l ), and 
r = v e + fi e k\ is the damping rate. Here ui is the wave frequency, k = kj_ = (k x ,k y ) is 
the wave vector, u = (c/B )d</>o/dx = — (c/B )E is the equilibrium electron drift in the 
presence of the dc electric field E , and k\ = k 2 + k 2 . In the absence of the dc electric field 
and dissipation, we observe from (3) that there exists frequency condensation at k y 3> k x . 
On the other hand, in a uniform plasma without Eq, we have the electron convective cell 
mode [3], to = —i{uj 2 e /uj 2 H ){v e + fi e kj_), which causes the cross-field electron transport due to 
random walk of the electrons in the electron convective cell electric field, even in a thermal 



equilibrium plasma 

Equation (2) describes the evolution of the electron fluid vorticity in 2D plane, e.g. 
in the (x, y) plane. The kinematic electron-thermal gyroviscosity (say due to turbulent 
fluctuations) has been introduced to accommodate the damping of fluctuations at short 
scales. 

To investigate the turbulence dynamics, we should have the knowledge of constant of 
motions of Eq. (2). The latter, without dissipation, admits two conserved quantities, namely 
the energy and enstrophy (squared vorticity) 

W = J(V ± (f)) 2 dxdy, Vl = J(V 2 ± (j)) 2 dxdy, (4) 

which show that the energy and mean squared vorticity are conserved ideally (or inviscidly) 
by 2D nonlinear interactions in a strongly magnetized electron plasma. 



III. SIMULATION RESULTS -FORMATION OF COHERENT STRUCTURES 

The modified NSE Eq. (j2J) has been integrated numerically with the help of a fully 
de-aliased pseudospectral scheme. Periodic boundary conditions are imposed along the 
x and y-directions. The electrostatic potential is discretized in a Fourier space using 
/(k, t) = J2it f{ r , t) exp(— ik ■ r). All fluctuations in our simulations are initialized with 
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FIG. 1: Simulation results from a random initial state lead to the formation of large scale coherent 
mode in potential fluctuations owing to an inverse cascade process (see left fig). Figure on the 
right shows the small scale eddies that are formed essentially from the forward cascades of electron 
fluid vorticity at a similar time. This is called a dual cascade phenomenon, observed ubiquitously 
in many 2D turbulent systems. Here we have chosen V /uj ce = 0.01 and (3 = {u)p e /{jj 2 H )K n L = 2, 
where L is the box size. 

a Gaussian random number generator to ensure that the Fourier components are all spa- 
tially uncorrelated and randomly phased. This ensures that the choice of initial state is 
highly isotropic, i.e. k x ~ k y at t = 0. Similarly, the boundary conditions (periodic in x, y 
directions) do not impose any kind of anisotropy. Moreover, the results to be presented here 
are independent of the size of computational domains, number of Fourier modes, as well as 
the integration time steps. 

The initially normalized energy spectrum, peaked at /c m i n , is chosen to lie within the 
wavenumber interval k m i n < k < k max /2. During the evolution of our simulations, the 
turbulence eventually evolves through vortex-merging in which like-signed smaller length- 
scale fluctuations merge to form relatively large-scale fluctuations. The process continues 
until all merging has occurred to finally form the largest scale coherent vortex dominated 
by the minimum allowed k in the simulation. The inertial range turbulent cascades in 
such a manner leads to the formation of large scale structures in potential fluctuations. By 
contrast, the vorticity structures break up and progressively form smaller eddies by means 
of a forward cascade processes. The large scale potential field nonetheless coexist with the 
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small scale vorticity fluctuations and lead essentially to a dual cascade process. This is 
shown in Fig. ([1]) when the turbulence has reached its saturated state. The formation of 
large scale potential and small scale vorticity structures through nonlinear interactions can 
be understood in the context of a dual cascade phenomenon as the 2D electron plasma 
system admits two inviscid invariants. Under this process the potential cascades towards 
longer length-scales, while the fluid vorticity transfers spectral power towards shorter length- 
scales. The randomly excited Fourier modes transfer the spectral energy in the inertial range 
by conserving the constants of motion (i.e. the inviscid invariants) in k space. This leads to 
a statistically stationary inertial regime associated with the forward and inverse cascades. 
The dual cascades in the electron plasma wave turbulence is similar to the one that occurs in 
2D hydrodynamic turbulence. The energy cascades towards smaller scales in the simulations 
is terminated essentially by a kinematic gyro-viscous damping. The latter is efficient at the 
smaller turbulent scales. The dual cascade processes reported here are further consistent 
with the number of inviscid constants (see Eq. (jH)) as admitted by the underlying system 
ofEq. (E]). 

Note that the presence of the density gradient convection caused by mean electric field 
does not quantitatively alter the characteristic nonlinear features of 2D electron plasma 
wave turbulence, as it exhibits a cascade property similar to that of the Charney-Hasegawa- 
Mima equation. The evolving relaxation of the turbulent fluid is independent of spatial and 
temporal resolutions, as well as higher turbulent Reynolds numbers. The latter slow the rate 
of relaxation, while the qualitative physics remains more or less unaltered. Interestingly, we 
notice the formation of thin current sheet like structures in the evolution of the vorticity as 
shown in Fig. (pQ). Furthermore, the formation of the observed large scale structures in our 
simulations can be understood from an exact nonlinear analytic solution of Eq. (T5]). This is 
described below. 

IV. ANALYTICAL SOLUTIONS FOR A DOUBLE VORTEX 

An analytic understanding of the emergence of coherent nonlinear structures observed 
in our fluid simulation can be achieved merely qualitatively by considering exact nonlinear 
solutions of the potential in Eq. ([2]). The underlying physical mechanism for the formation of 
coherent structures rests on the nonlinear interactions between different scalesize modes. In 
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saturated states, the vector product nonlinearity overwhelms the dissipation, and a balance 
between the wave dispersion and nonlinearity gives rise to long-lived stable vortex structures 
that we observed in our simulations. 

The two-dimensional traveling wave solutions of Eq. (j2j), without dissipation, in the 
moving frame can be obtained by supposing that x = x, £ = y — ut,t = t, where u being 
the speed of the vortex. This further transforms the fluctuating quantities into <p(x, y, t) = 
(p(x, £). The potential fluctuation obeys the localized boundary conditions — > as y — > oo 
for all x. The condition of localization associated with the coherent structures further enables 
us to use the periodic boundary conditions in computational domain of the nonlinear fluid 
simulations. The perturbed entities, or coherent structures, in the local region are far from 
the boundary and can translate freely. Hence no boundary effects are considered in our 
analysis. Boundary effects can nevertheless be important for certain processes, but this issue 
is well beyond the scope of the present work. Moreover, the x — £ plane in the simulations is 
orthogonal to the background magnetic field. On using the transformed co-ordinate systems, 
the electron vorticity Eq. (j2J) can be translated into a co-ordinate system that is moving 
along ^-direction with the vortex speed. Following the previous treatment of the vortex 
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dJvI<P--^—A=0, (5) 
V (u-uo) J 

where = {d/d£)-[c/{u-u Q )B Q }[{d(j)/dx)-{d(t)/d£)d/dx}, V\(f) = (d 2 (f)/dx 2 ) + (d 2 (f)/d^ 2 , 
and = —{ijJ 2 e uj ce /uj 2 H )K n . Equation (5) is satisfied by 

VI = C x <j> + C 2 x, (6) 

where C\ and C 2 are constants, and C\ + [(c/ (u — u )B ]C2 = u*/ (u — u ) = a 2 . For a 2 > 0, 
Eq. (6) admits a double vortex [l5]. The profiles of the latter in the outer (r > R, where 
r 2 = x 2 + £ 2 and R is the vortex radius) and inner (r < R; centered at the point x = £ = 0) 
regions are, respectively, 

(j)(r , 6) =(f) K 1 (ar) cos6, (7) 

and 

<p(r, 9) = \<Pi Ji ( 7 r) + Ar/7 2 1 cos 9, (8) 
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FIG. 2: Evolution of the total energy in the simulation. 

where cosfl = f/r, <j) = DiR/(a 2 + ^K^aR), & = -a 2 DiR/j 2 (a 2 + 7 2 )Ji(ai?), A = 
(a 2 + 7 2 )(w — uq)Bq/c, and the constant 7 has to be determined from the transcendental 
equation 

K 2 (aR) = HlR) , , 

aK x (aR) iM^R)' 

Here J\{J2) and K\{K2) are the Bessel and modified Bessel functions of the first (second) 
order, respectively. It turns out that the presence of the density gradient is essential for the 
formation of a double vortex. 

V. EVOLUTION OF TURBULENT ENERGY AND SPECTRUM 

The evolution of the turbulent energy associated with the potential and vorticity fields 
is shown in Fig. (j2j) where the linear and nonlinear phases of the evolution are clearly 
marked by the curve. The nonlinear Eq. (J2J), in the absence of dissipation, conserves 
the total energy. In the presence of dissipation, small scales are dissipated. Since the 2D 
electron plasma system is dominated by the smaller fc's, the large scale turbulent eddies 
tend to grow and contain most of the turbulent energy. The volume integrated energy 
therefore grows in time until the steady state is reached. The evolution of the energy is 
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FIG. 3: The 2D electron plasma turbulence in the presence of the density gradient (maintained by 
an external electric field) exhibits a turbulent spectrum close to the Kolmogorov-like k~ 5 ^ 3 scaling 
over the inertial range. The effect of a variation of the mean electric field is also shown in the 
turbulent spectrum. 

shown in Fig. (j2J). This is further consistent with the formation of large-scale coherent 
convective cells/flows in which the energy associated with turbulence must evolve in order 
to relax the turbulence into a well organized coherent structure. After nonlinear interactions 
are saturated, the energy in the turbulence does not grow and remains nearly unchanged 
throughout the simulations. Correspondingly, the energy transfer rate shows a significant 
growth during the linear and initial nonlinear phases (not shown in Fig. (j2J)). However, 
when nonlinear interactions saturate, the nonlinear transfer of the energy in the spectral 
space amongst various turbulent modes becomes inefficient and the energy transfer per unit 
time tends to become negligibly small. 

The power spectrum associated with the electron plasma turbulence exhibits a spectral 
slope close to a Kolmogorov-like scaling k~ 5 ^ 3 , as shown in Fig. ([3]). This is indicative of the 
eddy interaction being the dominant process in the spectral cascades of the inertial range 
turbulent energy despite the presence of the convection of the density gradient by the mean 
electric field, and is consistent with the Kolmogorov-like phenomenology. 
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FIG. 4: Evolution of the effective diffusion coefficient is studied in the presence of the density 
gradient in a nonuniform 2D electron magnetoplasma. Increasing strength of the inhomogeneity 
leads to a reduction in the net electron transport across the magnetic field lines. This result is also 
consistent with the experimental observations. 

VI. TURBULENCE TRANSPORT 



The cross-field electron transport is triggered primarily by the mutual collisions of the 
electrons in the presence of an ensemble of coherent structures. Predominantly, it is the 
successive collisions amongst the electrons that can potentially lead to a net outward drift 
across the external magnetic field. However, the presence of an external orthogonal electric 
field is believed to improve the stability of 2D electron plasma. Motivated by this observa- 
tion, we investigate here the transport of electrons that is influenced by turbulent evolution 
of the potential and vorticity in 2D electron plasma. We follow the evolution of the effective 
diffusion by increasing the magnitude of externally imposed mean electric field. An effective 
electron diffusion coefficient can be calculated from -D e // — /o°(V_L(r, i) ' Vj_(r, £ + t'))dt', 
where the angular bracket represents spatial averages. The perpendicular component of the 
electron fluid velocities in 2D electron plasma is V^ « z x V0. Since the 2D EDM tur- 
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bulence is confined in a plane orthogonal to the ambient magnetic field in z direction, the 
effective cross-field diffusion coefficient, D e ff, essentially describes the diffusion processes 
associated with the transverse motion of the electrons in x — y plane. We compute -D e // i n 
our simulations to measure the turbulent transport that is associated with the intermittent 
electron plasma wave turbulence. It is observed that the effective cross-field transport is 
lower, when the field perturbations are Gaussian. On the other hand, the cross-field diffu- 
sion increases rapidly with the eventual formation of longer length-scale structures. This is 
shown in Fig. (j3J), which exhibits a dependence of D e ff with time. Thus, an enhanced cross- 
field transport level results primarily due to an emergence of large-scale coherent structures 
in a nonuniform magnetized plasma. The enhanced cross-field diffusion coefficient observed 
in our simulations is, therefore, consistent with the generation of large-scale flows. 

The most remarkable point to emerge from our simulations is the reduction of the electron 
transport across the confining magnetic field when the strength of the applied or external 
electric field is increased. We have reported three different cases for the evolution of the 
diffusion coefficient corresponding to (3 = 0.2,0.5 and 2.0 in Fig. (jlj). It is clear from Fig. 
(J4j) that the effective transport associated with the electron plasma waves is much better 
(see the lower curve in Fig. (j3]) for (3 = 2.0) for a higher strength of the applied or external 
electric field. Our simulation results are further consistent with the experimental observation 
ofRef. J 



VII. DISCUSSION 



In this paper, we have studied the turbulent properties of nonlinearly interacting low- 
frequency EDMs in a nonuniform magnetoplasma by using computer simulations. Specif- 
ically, we have numerically investigated the modified Navier-Stokes equation (2) and have 
studied spectral cascades and cross-field electron transport produced by nonlinearly sat- 
urated turbulent modes in such a plasma system. It is found that self-organized coherent 
structures emerge due to nonlinear interactions between different scale sizes 2D fluctuations. 
The final state is thus intermittent and exhibits a dual cascade in which a forward cascade 
of the electron vorticity coexists with an inverse cascade of the potential fluctuation. The 
inertial power spectrum is close to the Kolmogorov-like 5/3 scaling, whereas the transport 
is dominated by large scale vortical structures. We find that the stability of the electron 
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plasma column confined by a straight magnetic field can be improved by virtue of an exter- 
nal electric field maintaining the density inhomogeneity, a result that is also consistent with 
experimental observations [2|. 
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